arXivil509.04283v2 [astro-ph.GA] 8 Feb 2016 


Draft version February 9, 2016 

Preprint typeset using DTeX style emulateapj v. 01/23/15 


THE PHOENIX STREAM: A COLD STREAM IN THE SOUTHERN HEMISPHERE 

E. Balbinot\ B. Yanny^, T. S. lP, B. Santiago^’^, J. L. Marshall^, D. A. Finley^, A. Pieres^’^, T. M. 

C. Abbott®, F. B. Abdalla^, S. Allam^, A. Benoit-Levy^, G. M. Bernstein®, E. Bertin^’^®, D. Brooks^, 

D. L. Burke^^’^^, A. Carnero Rosell®’^®, M. Carrasco Kind^^’^®, J. Carretero^®’^^, C. E. Cunha^\ L. N. da 
C osTA®’^®, D. L. DePoy®, S. DesaP®’^^, H. T. DiehU, P. Doel^, J. Estrada^, B. Flaugher^, J. Frieman^’^®, 

D. W. Gerdes^\ D. Gruen^^’^®, R. A. Gruendl®^’^®, K. Honscheid^^’^®, D. J. James®, K. Kuehn^®, N. Kuropatkin^, 
O. Lahav^, M. March®, P. MartinP^’^^, R. Miquel^®’^^, R. C. NichoU®, R. Ogando®’^®, A. K. Romer®®, E. Sanchez®\ 
M. SchubnelUP I. Sevilla-Noarbe®®’®^, R. C. Smith®, M. Soares-Santos^, F. Sobreira^’®, E. Suchyta^^’^®, 

G. Tarle^®, D. Thomas^®, D. Tucker^, A. R. Walker® 

(The DES Collaboration) 

Draft version February 9, 2016 


ABSTRACT 

We report the discovery of a stellar stream in the Dark Energy Survey (DES) Year 1 (YlAl) data. The 
discovery was made through simple color-magnitude filters and visual inspection of the YlAl data. 
We refer to this new object as the Phoenix stream, after its resident constellation. After subtraction 
of the background stellar population we detect a clear signal of a simple stellar population. By 
fitting the ridge line of the stream in color-magnitude space, we find that a stellar population with 
age r = II.5 ± 0.5 Gyr and [Fe/H] < 1.6 located I7.5T0.9 kpc from the Sun gives an adequate 
description of the stream stellar population. The stream is detected over an extension of 8?I (2.5 
kpc) and has a width of ^54 pc assuming a Gaussian profile, indicating that a globular cluster is 
a probable progenitor. There is no known globular cluster within 5 kpc compatible with being the 
progenitor of the stream, assuming that the stream traces its orbit. We examined overdensities along 
the stream, however no obvious counterpart bound stellar system is visible in the coadded images. 
We also find overdensities along the stream that appear to be symmetrically distributed - consistent 
with the epicyclic overdensity scenario for the formation of cold streams - as well as a misalignment 
between the Northern and Southern part of stream. Despite the close proximity we find no evidence 
that this stream and the h alo cluster NGC 1261 have a common accretion origin linked to the recently 
found EriPhe over density ( Li et al.|2QI5 ). 
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1. INTRODUCTION 

Our understanding of the structure of the Galactic halo 
has evolved considerably in the past two decades, largely 
thanks to deep and homogeneous photometric s urveys, 
such as the Tw o Micron All Sky Survey (2MASS; Skrut- 
skie et al.|2QQ6 ) and the Sloan Digital Sky Survey (SDSS; 
Ahn et al.|2Ul4|) . The stellar halo is now known to be in- 


habited by a variety of spatial and kinematic stellar sub¬ 
structure, from globular clusters (GCs) and dwarf galax¬ 
ies to extended stellar clouds and streams (see e.g. Will- 
man et ah 2005 Belokurov et ah 2006| 2007). In tact, 
recent simulations based on hierarchical models of struc¬ 
ture formation predict that most halo stars were brought 
by the disr uption of the Galactic substructures (Bullock 
etaLl[2M] ). 


I'he thin and cold stellar streams found in the Galaxy 
often span tens of degrees on the sky and originate from 
the tidal effects of the host on the progenitor, whether a 
GG or a dwarf galaxy. Perhaps the most conspicuous ex¬ 
amples of Galactic streams are those associated with the 
Pal 5 GC and the Sagittariu s dwarf (Odenkirchen et al. 
2001 Newberg et al. 2002). The tidal nature of such 
strearns makes them usefu l probes of the dark matter dis¬ 
tribution ac ross the halo (Johnston et al. 2005 Kiipper 
et al. 2015). Detailed modeling of a stream's position, 
distance, kinematics, gaps and overdensities in extended 
cold streams also leads to constraints on the amount of 
dark matter fragments orbiting the halo , known as sub¬ 
halos ( Yoon et al.|2011 Ngan et al.|2015 ), on th e progen¬ 
itor’s prop erties and on Galactic parameters (Koposov 

eraI][2MQ| ). __ 

I'he Dark Energy Survey (DES; |Abbott et al.||2005) 
is an on-going deep {g ^ 24.7) photometric survey m 
the Southern hemisphere that started its planned 5-year 
mission of collecting data in 2013. Despite its focus on 
cosmology, DES data have already produced a wealth 
of results pertaining to resolved stellar populations in 
the Galaxy and its vicinity, including the analysis of the 
structure and stellar po pulations in the outs kirts of the 
Large Magellanic Gloud (Balbinot et al.|2Q15|), the identi¬ 
fication of new Galac t ic companions ([Bechtol et al.|2015 ; 
Koposov et al.||2Q15l ILuque et al.|[2(jl5[ [Drlica-Wagn'^ 
et al. |2015| |K im & Jerjen|2015p, and the development of 
a new search for variable stars exclusively based on DES 
data (Hatt et ah, in preparation). 

We here report on the discovery of the first cold stellar 
stream using DES data. In Section 2 we give more details 
about DES, the data used, and the search algorithm. Our 
results are presented in Section 3 and our conclusions are 
in Section 4. 


2. DATA ANALYSIS 

DES is a wide-held optical imaging survey using broad 
photometric bands {grizY) performed with the Dark En¬ 
ergy Camer a (DECam; described in detail in Elaugher 


et al. 2015|). The DECam focal plane is coniprised of 
74 CCDs: 62 2kx4k CCDs dedicated to science imaging 
and 12 2kx2k CCDs for guiding, focus, and alignment. 
DECam is installed at the prime focus of the 4-meter 
Blanco telescope at Cerro Tololo Inter-American Obser¬ 
vatory. In this conhguration, DECam has a 2.2-degree¬ 
wide held-of-view and a central pixel scale of 0.263 arc- 
seconds. The full DES survey is scheduled for 525 nights 


distributed over hve years. Here, we consider data from 
the hrst year of DES obtained between 15 August 2013 
to 9 Eebruary 2014. 

The hrst internal annual release of DES data (YlAl) 
comprises the data products obtained from the process' 
ing of a subset of wide- and supernova- held d ata accumu 


lated during the hrst year of DES operations (Diehl et al 


2014). Briehy, the image processing pipeline consists of 


image detrending (crosstalk correction, bias subtraction, 
hat-helding, etc), astrometric calibration, nightly photo¬ 
metric calibration, global calibration, image coaddition, 
and object catalog creation. Eor a more detailed de¬ 
script ionof_the_^^D^£Jma£e_Qroce^ng_Qmeline, we 


refer to Desai et al. (|2012|); Mohr et ah] (12 012) and for 
a recent overview see |Balbinot et al.|(|2015p. The S Ex¬ 
tractor toolkit is used to create image catalogs from 


the processed an d coadded images (Bertin 2011 


Bertin 


fc Arnoiits||1996 ). The number of overlapping exposures 


in YlAl varies, but most of the footprint has at least 
3 coadded exposures. The YlAl coadd object cata¬ 
log contains ^ 131 million unique objects spread over 
^ 1,800 deg^. This area includes ^ 200 deg^ overlapping 
with the Stripe-82 region of SDSS, as well as a contigu¬ 
ous region of ^ 1,600 deg^ overlapping the South Pole 
Telescope (SPT) footprint. 

We perform stellar selection on the YlAl coadd ob¬ 
ject catalog based on the spreadjnod el quantity output 
of SExtractor (Desai et al. 2012). To avoid issues 
arising from htting the PS4' across variable depth coad¬ 
ded images, we utilize the weighted-average (wavg) of the 
spread-model measurements from the single-epoch expo¬ 
sures. Our stellar sample consists of well-measured ob¬ 
jects with \wavg spread-model-i\ < 0.003, flags < 
4, and magerr-auto-{g^r^i} < 1 (henceforth referred to 
as “stars”). Our stellar completeness is >90% down to 
magnit ude q ^ 22, at which point it drops to ^ 50% by 
g ^ 23 ( Bechtol et al.||2015 ) 


Stars are extinction corrected according to ISchlegel 
et al. J 19981) with t he scaling correction from jSchlany 


^ 44nkbeiner ( 2011|) assuming the extinction curve from 


Cardelli et al. ^1989 ) and a calibration at infinity, that 


is we assume that the light of every object in our sample 
crosses the full extent of the dust column measured in 
the dust maps. 

Erom this point all magnitudes considered in this paper 
are corrected for the extinction. 


2.1. Search method 

Using the objects classified as stars according to the 
criteria described in the previous section, we apply nar¬ 
row color filters to isolate interesting stellar types such as 
old turnoff stars and horizontal branch (HB) stars. To 
avoid issues related to the inhomogeneous photometric 
depth of the survey and saturation of bright stars, we 
only use stars with magnitudes 17 < ^ < 23. We find 
that this magnitude limit yields a sample that has a com¬ 
pleteness that is fairly constant across the footprint and 
produces a smoothly varying density for the field stars. 

Eor each color-selected catalog we build a sky “den¬ 
sity map”. Throughout this paper we use density maps 
to refer to maps where we show the number of sources 
per pixel (N) in a Cartesian projection. The pixel size is 
made explicit whenever necessary. The pixel area is cor- 








































































































Figure 1. YlAl density map for stars with 0.2 < (g — r) < 0.6 and 20 < < 23. The open triangles show the anchor points adopted for 

the stream. Other interesting objects are labeled in the figure. This density map was convolved with a 2 x 1-pixel Gaussian. Each pixel has 
a size of 4.5' x 2.8'. 


reeled for changes in solid angle with declination. These 
density m aps are visually inspected for over densities. In 
[Figure l| we show the particularly interesting density map 
for stars with 0.2 < {g — r) < 0.6 and 20 < < 23, 

which selects mainly turnoff and upper main-sequence 
stars from an old ( > 10 Gyr) simp l e steh ar population 


(SSP), according to Bressan et al. (2012). Several fea- 


tures are noticeable, su 
1261, the Phoenix dw; 
[I 977 I, and two of the : 

Lch as the globular cluster NGG 
arf galaxy (Ganterna & Elower 

recently discovered dwarf galax- 

ies ([Bechtol et al.| 

TT U? 1 1 \ 0 /H Va 

2015 

Koposov et al. 

TO II 1 I 1 \ 

|2015 

): Reticulum 


labeled with their names in the figure. A linear struc¬ 
ture is also visible near the Phoenix dwarf extending 
from {a, 6) - (20°,-57°) to (27°,-45°). This struc¬ 
ture is highlighted by open triangles marking high den¬ 
sity points along the stream candidate. 

In the same figure, a large overdensity of stars is visible 
between the stream candidate and NGC 1261. This fea¬ 
ture is the Eridanus-Phoenix (EriPhe) overdensity and 
it is discuss ed in detail in a simultaneous publication (O 


3. RESULTS 

By means of the method outlined in the previous sec¬ 
tion we perform a visual search for stellar overdensi¬ 
ties. This search, conducted over the full YlAl footprint, 
ha s revealed only on^ stream candidate which is shown 
in [Eigure l] Eor simplicity, we refer to this candidate 
stream as the Phoenix stream due to its proximity to 
the Phoenix constellation. 

To study the stellar population that comprises the 
Phoenix stream we define a line passing through the cen¬ 
ter of the stream using anchor p oints along the stream 
(three open triangles in Eigure 1). We then select stars 


^ Recently [Mackey et al.|(|2015|) reported a 10 kpc long stream 
associated to tJie Large Magellanic Cloud (LMC) using a similar 
data set to the one in this work. We confirm this detection, however 
its tidal origin is still uncertain. 


inside a box defined by the stream central line with an 
offset of ± 0.8° in RA. We name this selection on stream. 
To compare with the typical Milky Way (MW) stellar 
population at this position in the sky we select stars in 
boxes that are offset by ±1.5° with respect to the central 
one and have the same width as the on stream region. 
We name these selections off stream east and west. Eor 
each region described above we compute the solid an- 
gle normalized Hess diagram. We use Mangle masks 
( Swanson et al.|[2QQ8 ) to compute the solid angle of each 
box takin g into acc ount possible holes in the survey foot¬ 


print. In Eigure 2] we show the Hess diagram of the on 


stream minus the average diagram of the two off stream 
ones in logarithmic scale. _ 

Erom the decontaminated Hess diagram shown in [Eig-[ 
lure 21 we estimate that the stream has ^ 500 stars that 
fall within the photometric limits of DES. This decon¬ 
taminated Hess diagram was smoothed using a Gaussian 
kernel with a dispersion of 0.06x0.2 in color and magni¬ 
tude respectively. This step is required due to the low 
number of stars in the stream and allo ws ns to define 
a ridge line shown as the red circles in [Eigure 2] The 
ridge line is defined as the peak value of counts in color 
for each magnitude bin, the bin size for the ridge line 
construction being twice as large as the one used for the 
Hess diagram. Magnitude bins with low counts or peak 
values that obviously depart from the bulk of the stream 
stars are discarded. We define the error bars as equal 
to the Gaussian kernel size, which is always larger than 
the photometric errors in the magnitude range shown. 
This choice of error bar accounts for the broadening of 
the Hess diagram due to the smoothing process. Using 
the typical photometric error would yield unrealistically 
small uncertainty estimates that would propagate into 
the stellar evolution model fitting, described below. 

We compare the ri dge line to different PARSEC stel¬ 
lar evolution models (Bressan et al.[[2012) by computing 
the minimum distance from each ridge line point to a 
given model. The model grid has a resolution of 0.01 in 





































described by a model with {m — M) = 16.21 ±0.11 (or 
dQ = 17.5 ± 0.9 kpc), logio(r/^r) = 10.06 ± 0.02 (or 
r = 11.5±0.5 Gyr), and Z < 0.0004 (or [Fe/H] < -1.6). 
The best-fit model is shown in [Figure'^ as the solid black 
line. Notice that the lowest metallicity available in our 
model grid is still consistent with the stream CMD, thus 
we are only able to define a upper bound for the metal- 


licity. We summarize the stream parameters in Table 1 
We use the best-fit model to define a region in the color 
magnitude diagram (CMD) where stream stars are more 
likely to be. This region is shown as the dashed line in 


[Figure 2 The locus shown in the figure was defined by 
color shifting the best-fit model by twice the typical color 
error at each magnitude value. We also consider that the 
color uncertainty at magnitudes brighter than ^ = 21 is 
cons tant and e qual to 0.03. 

we show the density map for stars (left 


Figure 3 

) built using the color selection described above but 


Figure 2. Decontaminated Hess diagram of the stream candidate. 
The decontamination process is described in detail in the text. The 
solid line shows an 1 1.5 Gyr and [Fe/H]=-1.9 PARSEC isochrone 
(jBressan et al.|2012|). The dashed lines show the color magnitude 
diagram (UlViDj region selected to isolate stream stars, including 
a box to select Horizontal Branch (HB) stars. The red dots with 
error bars form the ridge line used to perform the isochrone fit. 


\ogiQ{T/yr) in the range from 9 to 10.16 and 0.0002 in 
Z in the range from 0.0001 to 0.001, where r and Z are 
age and metallicity. The distance modulus was explored 
in the range from 15 to 18 in steps of 0.01. For each 


In 

pane' 

only for stars with 20 < ^ < 23. The best-fit model de¬ 
scribes the stream population in the full domain of colors 
and magnitudes observed, however we find that the HB, 
Sub-giant Branch (SGB), and Red Giant Branch (RBG) 
are very sparsely populated. Including these stars in the 
GMD selection adds more noise than signal to our density 
maps. In the same figure, we also show the density map, 
using the same selection as before, but for sources clas¬ 
sified as galaxies (center panel ). On the rightmost p anel 
we show a reddening map from Schlegel et al. (1998). On 
the last two panels the solid line shows the position of 
the stream. We notice no obvious features in the galaxy 
or reddening distribution that could mimic the presence 
of the stream. 

3.1. Possible progenitors 

In order to investigate any possible progenitors for the 
new stream we assume that str eams ye ap proximate 
tracers of the progenitor’s orbit (Bovy 2014). We also 


exploit the property of spherically symmetric potentials 
in which orbits should be confined to a plane contain¬ 


ing t he center of such a potential (Binney & Tremaine 


2008) 


parameter combination we compute the probability that 

ric potentials ( 

Johnston et al. 

1996| 

a given ridge line point was drawn from the isochrone at 

potentials this assumption onJ 

[y hol( 


its minimum distance position to that given point. The 
probability is computed assuming a normal distribution 
with a Icr dispersion as indicated by the error bars. All 
ridge line points are given the same weight. An obvious 
improvement would be to weight these points by a mass 
function (MF). However, the MF of a stream is likely very 
different from an initial MF (IMF) and can vary along 
the stream itself ( Koch et al.||2004 ). For this reason, we 


For more complex 
y holds close to the pro¬ 
genitor itself. There is evidence that the MW potential, 
at least in its inner parts, i s well approximated by an 


axisymmetric potential (see Kiipper et al. 2015). The 


Phoenix stream lies at 18.4 kpc from the Galactic center 
where the MW potential should be reasonably spherical 


leave the study of the IVH' of the stream to future works 
with deeper photometry and more accurate membership 
probabilities. 

We calculate a likelihood function for our model by 
multiplying the individual probabilities of each ridge line 
point. We define the best model as the one maximizing 
the likelihood and the parameters’ uncertainties are de¬ 
rived using the profile likelihood technique (e.g. Rolke 
et al.||2005|). To estimate the 90% confidence interval of 


each htter parameter we find the value of that parame¬ 
ter where the log-likelihood (maximized with respect to 
the other parameters) decreases by 2.71/2 from its max¬ 
imum value. We find that the stream population is well 


(Kiipper et al.|2bl5t|Deason et al.|2011t|Bell et ai.|2008|). 
U nder the assumption that the strearn formed through 
the interaction with the MW potential only, we expect 
that it should be confined to a plane passing through the 
center of the Galaxy. When observed from the center of 
the MW this plane is described by a great circle. 

To define such a plane we choose three anchor points 
along the stream. These points are defined by their 
Galactic coordinates and the heliocentric distance to the 
stream. We find the plane that contains the anchor 
points and the MW center. And finally we find the circle 
oriented the same way as the plane that intersects all an¬ 
chor points. In order to intersect all three anchor points 
we must apply corrections to their heliocentric distances, 
which were so far considered all identical. 

We find that a heliocentric distance gradient of ~ 1 
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Figure 3. Left panel: density map after applying the color-magnitude selection shown in Figure 2. The white arrow points towards the 
MW center, eenter panel: density m ap using the same fil ter a s the previous panel but buil t using sources classified as galaxies. Right 
panel: E(B-V) color excess map fromlSchlegel et al.|(|1998|) andlSchlafly & Finkbeiner| (|2011|). On the last two maps the solid line shows 
the position of the stream as seen in the stellar density map. All maps use a pixel size ot 4.5'^ X 2.8'. 


kpc is necessary for a circular orbit to intersect all three 
anchor points. In order to check this possible distance 
gradient we build two separate decontaminated Hess dia- 
grams following the same procedure outlined in |section~3 


one using stars north of (5 = —56° and the other using 
stars South of this same declination value. Using the 
ridge line of each of these new Hess diagrams we pro¬ 
ceed with the fitting process, however this time we keep 
the metallicity and age fixed at the best-fit values found 
previously using the full stream length. We find that 
the best-fit distance modulus for the North part of the 
stream is 16.19 ± 0.12, and 16.35 ± 0.12 in the South. 
This is consistent with, but does not require, the ^1 kpc 
gradient required for a circular orbit. We conclude that 
a distance gradient cannot be ruled out for the stream. 
Detailed spectroscopic observations must be used to iso¬ 
late stream member stars based on radial velocities and 
chemical composition and confirm this scenario. For the 
purpose of looking for possible progenitors close to the 
stream we will assume a circular orbit. 

Using the best-fit circular orbit described above, we 
look for possible known GCs that are not in the YlAl 
footprint that could be progenitors. This approach does 
not explore other kinds of orbits, which are more likely, 
however it provides a useful approach to se arch for p ro- 
genitors in the close vicinity of the stream. In |Figure~^ we 
show the best-fit great circle in an all-sky Ait oft’ Galactic 
projection. We also show known GCs with Galactocen- 
tric distances between 15 and 25 kpc, distances consistent 
with the stream distance. We find that no GC is consis¬ 
tent with this stream under the assumption of a circular 
orbit. 

Another possible scenario is that the progenitor has 
completely dissolved and only its remains are visible 
along the stream. To investigate possible progenitors for 
this stream we first look for overdensities on the stream. 
We start by creating a reference frame where the hori¬ 
zontal axis is orient ed along the stream, si milar to what 
has been adopted by Majewski et al. (2003) for the Sagit¬ 
tarius stream. To create such a reference frame we use 
two Euler angles that define two consecutive rotations 
(0, 6>). The angles were determined by finding the plane 



Figure 4. Aitoff projection in Galactic coordinates. The gray 
circles show GCs with Galactocentric distances between 15 and 25 
kpc. The solid line shows the great circle that best fits the Phoenix 
stream, the portion of the stream observed is highlighted with a 
broader line. 


that intersects the anchor stream points in Equatorial 
coordinates. The new reference frame has an azimuthal 
component (A) with an arbitrary origin and defined in 
the range [0, 27r), and an elevation component (/3) defined 
in the range [— f, f ]• The values of cj) and 0 adopted are 
-29?698 and 72?247 respectively. 

In |Eignre~^ we show the density map of color-selected 
stars in the coordinate system described above. We also 
show the average density map in projection onto the P 
{left panel) and A {top panel) axes. The A projection is 
shown for stars that he within a = 0°.18 from the stream 
centroid in /d, where a is the standard deviation with 
respect to the stream median line. The stream median 
line is defined in the stream coordinate system as the 
peak in the P projected density. We also show the P 
projection for stars south (green dashed line) and north 
(blue dashed line) of the central overdensity. We notice 
an offset of ^ 0°.14 in P between the North and South 
portion of the stream. 

The stream width depends on two sets of factors. The 
first and more obvious is the progenitor’s size and veloc¬ 
ity dispersion. The second is the shape of the gravita¬ 
tional potential. Eor instance, triaxial potentials t end to 


increase the fanning of stream stars sig nificantly (Price- 
Whelan et al.|2Q16 Pearson et al.|2Q15 ). However, there 


is evidence that the inner halo of the Galaxy is relatively 
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Figure 5. Bottom right panel: density map in the stream’s coo rdinate sy stem built using color-magnitude selected stars. The white arrow 
points towards the MW center. The pixel size is the same as in [Figure l| The dashed lines shown in Figure 2 encompass the region of the 
CMD where stars were selected from. The left panel shows the average density map value in projection onto the (3 axis, the solid horizontal 
line marks the position of the stream centroid and the dotted lines shows the Tier limits. The green (blue) dashed line shows the same 
kind of /3 projection but for the South (North) part of the stream. The top panel shows the average density map value in projection along 
the A axis. This latter is built using only stars that are within Tier of the centroid of the stream. Vertical solid lines mark the position 
of candidate progenitors. In this panel the dashed line shows the typical background contribution, computed using stars within Tier of a 
center line offset 5cr from the stream’s original center line. 


Table 1 

Phoenix stream parameter summary 


Name 

Value 

Unit 

Description 

r 

11.5 ±0.5 

Gyr 

Age 

[Fe/H] 

< -1.6 

dex 

Metallicity 

dQ 

17.5 ±0.9 

kpc 

Heliocentric distance 

doc 

18.4 ±0.9 

kpc 

Galactocentric distance* 

G 

54 

pc 

Stream width^ 

( q ;, S)start 

(20, -57) 

deg 

Stream begin point 

(q:, (5)en(i 

(27,-45) 

deg 

Stream end point 

© 

8.1 

deg 

Stream length 


* Computed using Rq = 8.3 kpc; f assuming a Gaussian profile. 


spherical dKiipper et aD 2015|) out to ^ 20 kpc, thus 
th< 


allowing us to assume the stream width maps only to 
the progenitor size and velocity dispersion. Using the 
stream coordinate system we determine that it has an 
on-sky width of cr = 0°.18, which translates to ^ 54 pc 
at its distance. Typically, 50 pc is consistent with the 
tidal radius of the MW halo GCs. The fact that the 
stream forms a thin coherent structure several kpc long 
makes it plausible that the progenitor was in fact a GC. 

From the density projected in the A coordinate we ex¬ 
plore the presence of overdensities (ODs) as pos sible pro- 
genitors. We label these ODs in the top panel of [Figure"^ 
First we call attention to G1 and C2, which have a slight 
offset with respect to each other in the (3 direction, show¬ 
ing hints of trailing and leading tail misalignment with 
respect to the orbit (|Bovy|2014). Apart from G1 and C2 
we find two other peaks, one to the north (N) and another 


to the south (S) of the central ODs. These ODs stand 
out when com pared to t he typical background counts 
(dashed line in |Figure 5|). We compute the typical lo¬ 
cal background noise at the position of the north (south) 
OD by taking the standard deviation of the background 
counts in the north (south) portion of the stream. We 
find that both ODs (N & S) peak densities stand out 
more than 4cr with respect to the background. The sig¬ 
nificance values are listed on Table 2. 

The fact that ODs N and S are approximately equally 
separated from the central OD could poi nt to epicyclic 
Qverde nsities such as t he ones reported in |Kupper et al.| 
(2015) for Palomar 5. Table 2 summarizes the positions 
and angular separation of the overdensities with regard 
to the central peaks. From this table we observe that 
ODs to the north are systematically at higher /3 than 
those in the South, hinting to the misalignment men¬ 
tioned above. 

Using the misalignment of the Northern and Southern 
portions of the stream (hinted injFigure 5), we may infer 
from geometrical considerations alone that its Northern 
part is closer to the MW center, hence being formed by 
stars that leave the progenitor through the inner part of 
its orbit, forming a leading tail. By construction, the 
South portion forms the trailing tail. From this argu¬ 
ment, we conclude that the stream is moving from south 
to north. 

All overdensities were visually inspected on the coad¬ 
ded images and catalog, however we could not identify 
any stellar system (e.g. globular cluster) that might have 
given origin to the stream. This result is very puzzling. 



































Table 2 

Overdensity positions 


OD 

RA 

deg 

Dec 

deg 

1 

deg 

b 

deg 

A 

deg 

/3 

deg 

AC1C2 

deg 

a 

Cl 

23.75 

-49.89 

285.68 

-65.76 

287.50 

0.00 

- 

5.1 

C2 

23.77 

-50.40 

286.16 

-65.28 

286.94 

0.22 

- 

5.0 

N 

26.15 

-46.34 

277.74 

-68.10 

283.16 

-0.09 

4.10 

7.1 

S 

21.23 

-53.70 

292.10 

-62.72 

291.11 

0.03 

4.17 

5.0 


especially if the scenario described above is to be con¬ 
firmed. The fact that no progenitor is found, but classic 
signs of cold tail formation are observed could indicate 
that a progenitor was fully disrupted very recently. 


4. CONCLUSIONS 

We report the discovery of a stellar stream in the 
Southern hemisphere. Through the visual fit of stellar 
evolution models we found that this stream is comprised 
of an old (11.5T0.5 Gyr) metal-poor ([Fe/H] < —1.6) 
population that is 17.5 kpc away from the Sun and 18.4 
kpc from the Galactic center. Though close in projec¬ 
tion, the Phoenix stream is not related to the Phoenix 


galaxy, whi ch lies at 440 kpc from the Sun (Karachentsev 
et al.|2004D. 

I'hrough the extrapolation of the stream outside the 
YlAl footprint, we found no known GC that could be 
its progenitor; however more eccentric and/or non-planar 
orbits were not considered. 

We also investigate the distribution of overdensities 
along the stream in search for a progenitor. We find 
that none of the CDs has any obvious stellar overdensity 
associated with it when coadded images were inspected. 
We find that the CDs with high significance display a 
symmetric pattern with respect to a central overdensity. 
This central overdensity shows some hints of misalign¬ 
ment perpendicular to the orbit direction, which could 
indicate the position of the progenitor. 

A diffuse stellar overdensity has recently been found 
in the D ES data which nearly overlaps with the Phoenix 
stream (Li et al. 2015) . This overdensity (EriPhe ) was 
previously hinted at by Carballo-Bello et al. (|2014) as an 
anomalous background population close to NGG 1261. 
EriPhe and NGC 1261 share a similar heli ocentric dis- 
tance as the Phoenix stream. Using galpy (Bo vy||2Q15|) 
and l iterature proper motions for NGC 1261 (jUambis 
2006), we integrate the cluster orbit and find that it 


roughly aligns with the Phoenix stream and that its mo¬ 
tion is retrograde with respect to the Solar motion. The 
close proximity of NGC 1261 and the orbit alignment 
with the stream may suggest that they could share a com¬ 
mon origin with the EriPhe overdensity. However, the 
stream appears not to be in a retrograde orbit, favouring 
a scenario where the stream is independent from EriPhe 
or NGC 1261. Eor an extended discussion about this 


scenario we refer to Li et al. (2015). 
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